Universal response of quantum systems with chaotic dynamics 
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The prediction of the response of a closed system to external perturbations is one of the central problems in 
quantum mechanics, and in this respect, the local density of states (LDOS) provides a deep description of such 
a response. The LDOS is the distribution of the overlaps squared connecting the set of eigenfunctions with the 
perturbed one. Here, we show that in the case of closed systems with classically chaotic dynamics, the LDOS is a 
Breit-Wigner distribution under very general perturbations of arbitrary high intensity. Consequently, we derive 
a semiclassical expression for the width of the LDOS which is shown to be very accurate for paradigmatic 
systems of quantum chaos. This work demonstrates the universal response of quantum systems with classically 
chaotic dynamics. 
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The action of a perturbation on eigenfrequencies and 
eigenfunctions of a quantum system has been a subject of 
paramount importance since the beginning of quantum the- 
ory. Its understanding is at the heart of fundamental problems 
of quantum mechanics like dissipation, phase transition or ir- 
reversibility. The usual perturbation theory is a good start- 
ing point to describes successfully this effect when the per- 
turbation is small. However, approximated theories usually 
fail for strong perturbations and highly demanding computa- 
tional methods are needed to describe characteristics of the 
perturbed system. 

The local density of states (LDOS) or strenght function 
is a widely studied magnitude to characterize the effect of 
perturbations on quantum systems and has been extensively 
computed for different systems and perturbations |Q]-|3]. The 
LDOS is a distribution of the overlaps squared between the 
unperturbed and perturbed eigenstates. Let us consider a 
chaotic one parameter dependent Hamiltonian H(x), and its 
quantum counterpart H(x) with eigenfrequencies 0Bt(x) and 
eigenstates \j(x)). Then, the LDOS of an eigenstate |i(jto)) 
(that we call unperturbed) is given by, 

Pi(o),8x) = £ | (j(x) \i(x ))\ 2 S{co - coij), 

i 

with 8x =x — Xq and flfly = C0j(x) — <X>i(xo). Furthermore, to 
avoid a dependence on the particular characteristics of the 
state |i(jto)}, an average over n unperturbed states in a small 
frequency window is performed 



p((a,8x) = -Ypi{co,8x). 



(1) 



The LDOS is related with important measures of irre- 
versibility and sensitivity to perturbations in quantum systems 
as the survival probability and the Loschmidt echo (LE) Jfjl- 

. In fact, the LDOS is the Fourier tranform of the survival 
probability [9] and its width gives the decay rate of the LE 
for a small enough strength of the perturbation 0, 0- In 
this letter these relations are exploited to show that LDOS 
has Lorentzian shape, usually called the Breit-Wigner distri- 
bution, under very general perturbations of arbitrary high in- 
tensity. Moreover, we derive a semiclassical expression for 



the width of the LDOS, G sc , in chaotic systems. The de- 
rived expression only depends on the perturbation, while the 
properties of the system are taken into account through a uni- 
form measure in phase space. Although O sc is derived for 
local perturbations we show that it also works in the case of 
global perturbations as a consequence of the Lorentzian char- 
acter of the LDOS, and the requirement of statistical inde- 
pendence between perturbed and unperturbed eigenfunction 
sets. Of course, such a requirement imposes restrictions on 
the admitted perturbations and we discuss this point at the end 
of the letter. We test the ability of (7j C to predict the width 
of the LDOS in perturbed cat maps and the Bunimovich sta- 
dium billiard with boundary deformations, observing that it 
works very well in both systems, even for strong perturba- 
tions far away from the Fermi Golden Rule regime. These 
results demonstrate for the first time the universal response of 
quantum chaotic systems to perturbations of classical nature. 
The Fourier transform of Eq. (fl} is given by, 

&\p\{t,8x) = ~Y,e- imMt (i(xo)\e^ R \i(x )), (2) 
n 

where the sum runs over the amplitude fidelity of eigen- 
states (whose square modulus is the survival probability). Let 
us evaluate the previous sum semiclassically. Vanlcek has 
proposed an approximation of the amplitude fidelity, named 
the dephasing representation flUl . by assuming a classically 
small perturbation in such a way that the shadowing theorem 
11 1 Ofl is valid. Using such an approximation, Eq. (O leads to, 



&\p](t,8x)f*s J dqdpW{q,p) exp[-iAS, (q,p, Sx) /«] , (3) 

where AS t (q, p,8x) is the action difference evaluated along 
the unperturbed orbit starting at (q,p) that evolves a time t. 
Moreover, W(q,p) — (l/n)£W,(<7,/?), with being the 

Wigner function of \i(xo)). In chaotic systems, W(q,p) re- 
duces to a uniform distribution. 

In the case of local perturbations, the right hand side of Eq. 
OJ has been evaluated on a Poincare surface of section by 
Goussev et al. Ill 211 resulting, 



&\p](t,8x)fxe 



-w\ 



(4) 
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FIG. 1: (Color online) Width cr of the LDOS as a function of the 
scaled perturbation strength % = NSk, for a local perturbation. Solid 
symbols correspond to the quantum case and solid lines to the semi- 
classical calculation . We use N = 800 and go = 0.01. In panel 

(a) we use a width /3 = 0.2 (circles) and j8 = 0.4 (squares). In panel 

(b) p = 0.7. Left inset: Schematic figure showing the used local per- 
turbation. The scaled shear [f(q) = 27ze(q,k) /k, (with A: the strength 
of the perturbation)] is plotted as a function of q. The limits of the 
perturbed region are indicated with qo and q\ , being j8 = q\ — q$ its 
width. Right inset: The LDOS p for j8 = 0.7, with % = 78.4 and 
89.6. In solid (red) lines we plot periodized Lorentzian functions 
with the corresponding widths. 
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The average is evaluated on the region of surface of section 
where the local perturbation acts, and AS(q,p, 8x) is the ac- 
tion difference after one step on the surface of section. For 
numerical comparison we consider a rectangular region, 



,-iAS(q,p,Sx)/h\ = 1 



1 fP2 f12 



I 



a 



-iAS(q,p,8x)/n 



dqdp, (6) 



where p\, p2, q\ and qi are the limits of the perturbed region 
with area a = (p2 — Pi )(<7? — ffi )■ Moreover, rj, named the 
classical decay rate in Ref. 1 121 . is the probability to reach the 
perturbed region per unit time, 



a 



(7) 



where srf is the area of the Poincare surface of section and T 
is the mean mapping time. 



The inverse Fourier tranform of Eq. (0]) is a Lorentzian 
function, the so-called Breit-Wigner distribution, 



P(g),8x)&L(y,g>) = 



7t(a> 2 



(8) 



and we define its width as half the distance around the mean 
value that contains the 70% of the probability (actually, this 
value is only relevant for the numerical computation). Then, 
the semiclassical approximation of the width results, 



<7 cr = tan 



1.9637. 



(9) 



We stress out that the semiclassical approximation was de- 
rived in the limit of a — > 0. However, our final expression 
(Eq. (0) can be extended to arbitrary values of a as a re- 
sult of the Lorentzian character of the LDOS and the property 
of short range correlation of chaotic eigenfunctions. To clar- 
ify this point, let us consider the following basis sets: {I/' ')} 
be the set of unperturbed eigenfunctions, {I/' 1 ')} be the set 
resulting after applying a local perturbation over an infinitesi- 
mal region 8(X\, and {I/' 2 ')} be the set resulting after applying 
the perturbation 8(X2 to the previous system (the one with the 
perturbation 8(Xi). For 8a,\ and 8(X2 being perturbations over 
disjointed regions of phase space, their corresponding LDOS 
should be statistical independent because chaotic eigenfunc- 
tions have correlations of short range in phase space. There- 
fore by assuming independence, the LDOS connecting the 
first and third basis sets is simply derived by convoluting the 
previous ones, 

/ L(8y x , (o x )L(8y2, (0 - (Oi)d(Oi = L(8yi + 8y 2 ,(o). 

The new LDOS is also a Lorentzian function, with 8y\ {8j2) 
being the parameter corresponding to the first (second) pertur- 
bation. Then, by using Eqs. (0, © and (|7]i it is easy to see 
that 5/i + 8j2 1S just the parameter resulting from the pertur- 
bation 8(X\ + 8(X2- Following this procedure, we can now add 
a perturbation 8aj and so on up to fill a finite region a. We 
would like to stress that even though the property of statistical 
independence is very reasonable, there are exceptional pertur- 
bations where such a property is not satisfied; see the end of 
the letter for a discussion of this point. 

Let us show the power of the semiclassical approximation 
to describe the LDOS in quantum maps, where the mapping 
time is fixed to T = 1, the phase space is a torus of area 
= 1, and the Hilbert space has finite dimension (with 
2nh = 1 /N). We consider the cat map, a canonical example 
in classical and quantum chaos studies, perturbed with a non- 
linear shear in momentum, 
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■e(q,k) 



(mod 1), 



where e(q,k) = (k/2n)[cos(2Kq) — cos(4nq)}, with k be- 
ing the strength of the perturbation. The action differ- 
ence for one iteration of the map is given by AS(q,8k) — 
(8k/4n 2 )[sm(27tci) -sinf4^q)/2] |fH. 
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FIG. 2: Idem Fig.Q]for global perturbations. The perturbation in the 
main plot is the same as in Fig.fj] For the inset, the perturbation is a 
shear in momentum and position (see text for details). 

For local perturbations lfl4h . the shear £(q,k) is only ap- 
plied to a q window from qo to q\, with width j3 = q\ — qo 
[see left inset in Fig.Q"); so, a = j3. Furthermore, we take into 
account the fact that the spectrum of the cat map is periodic 
because of a compact phase space. This periodicity changes 
the form of the LDOS because the probability that leaves from 
one border returns to the other. By assuming no correlation 
between the existing and returning probabilities, the LDOS 
transforms into a periodized Lorentzian function 

L(%,o))= £ L(r,a-2jcj/x), 

7=-°° 

with cot being the variable that specify the spectrum of eigen- 
phases. This distribution provides a new relation between the 
width of the distribution and y, whose first correction with re- 
spect to Eq. (0 is Gse ' — Osc[l - (Osc/tc) 2 )]. On the other 
hand, a numerical computation reveals a linear term which 
should be related to correlations between the existing and re- 
turning probabilities. We obtain the following estimate, 

oi p] ra cv[l + 0.24cv - {a sc /nf\ , (10) 

by fitting the linear term to the numerical data. 

FigQ] shows the width, a, of the LDOS and its semiclas- 

sical approximation, trie , for three different windows in po- 
sitions where the perturbation is applied. We plot the width 
as a function of the scaled perturbation strength, % = N8k, in 
such a way that figures are insensitive to N. As it can be seen 
the semiclassical approximation works very well. Moreover, 
the right inset of FigQ] (b) shows that the LDOS is a peri- 
odized Lorentzian function even for strong perturbations far 
away from the Fermi Golden Rule regime, which is identified 
with the quadratic behavior close to the origin. 

Fig. |2] compares the width of the LDOS with aJe m the 
case of global perturbations (a = 1). The perturbation in the 
main panel is the one used previously while in the inset, we 
use the same shear in momentum plus the following shear 
in position: e(p,k) = — (k/2n)[sm(6l[p)/3 + cos(4ji/>)/2]. 
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FIG. 3: (T as a function of the perturbation strenght Sx for the stadium 
billiard in the region of wave number around 200 (solid circles). The 
semiclassical approximation a sc using Eq. (9} is plotted with solid 
line. Inset: two Bunimovich stadium billiards with different shape 
parameter; the perturbed system is indicated in dashed line. 



In the latter the action difference is, AS(q,p,8k) — 
(8k/4n 2 )[sm(27iq) - sm(4nq)/2 + cos(67t/>) - sm(4np)/2\. 
Differences near the peaks at % w 20 and % rj 50 are related to 
the poor accuracy of Eq. ( TTOb close to the saturation, which is 
given by the width of the uniform distribution, 0.7n « 2.20. 

To further demonstrate the power of the proposed semi- 
classical approximation in a realistic system, we consider the 
desymmetrized Bunimovich stadium billiard with radius r and 
straight line of lenght a. This system is fully chaotic lfl5ll 
and has great theoretical and experimental relevance. The bil- 
liard is perturbed by the boundary deformation displayed in 
the inset of Fig [3] The area of the billiard is fixed to the value 
1 + it/ 4, so the boundary only depends on the shape parame- 
ter x = a/r. The boundary deformations are parametrized by 
r(q,8x) — ro(q) +z(q,8x) n, where q is a coordinate along 
the unperturbed boundary ro(q) defines ^ and n is the out- 
ward normal unit vector to 'rf at rpiq) (an explicit expression 
for z(q, 8x) is provided in Ref. OlaO . We consider the usual 
Birkhoff coordinates to describe the classical dynamics of the 
particle; that is, the variables q and p = sin(0), with 6 the im- 
pinging angle with n. To compute a sc , the action difference 
between the unperturbed and perturbed orbit results in ITLzTl . 

AS(q,9,8x) = |p|AL = 2|p|z(4,<5x)cos(0), 

where AL is the lenght difference between the unperturbed 
and perturbed orbits, and p is the momentum of the particle. 
The mean time between bounces with the boundary is given 
by T = mwA/(|p| IU7I1 . with m the mass of the particle, 
A the area of the billiard and 3? its perimeter. Then, the de- 
cay rate results i n 77 = |p| j3/ (mnA), with j3 the width of the 
perturbed region llOl : for the selected perturbation, /3 = . 
Fig. [3]compares the numerical and semiclassical calculations; 
h = 1 and m — 1/2 are used. The width computed with the 
exact eigenstates is plotted with full circles and the semiclas- 
sical aproximation, a sc , is plotted in full line. The calculations 
displayed in Fig. [3]were performed around the wave number 
200, where the eigenstates of the billiard were computed using 
the scaling method |Q . The agreement between the quantum 
and the semiclassical calculation is excellent. We notice that 
while the full quantum computation of a in Fig. [3]is very time 
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consuming [r»7x 10 7 seg in an CPU Intel Core 2 6400], the 
semiclassical calculation is a simple two variable integral. 

One final point is to discuss the character of the perturba- 
tion in order to satisfy the used property of statistical indepen- 
dence between the perturbed and unperturbed eigenfunction 
sets. Let us consider as an extreme example a chaotic Hamil- 
tonian of the form kinetic plus potential energy, and where 
the perturbed Hamiltonian is obtained by a displacement of 
the potential energy. In this situation, the two systems have 
the same spectrum and the corresponding eigenfunctions are 
connected by the displacement. So, the two sets are strongly 
correlated and the LDOS does not satisfy the Breit-Wigner 
distribution; see for instance Ref. B19I1 where this type of per- 
turbations is analyzed. From the classical point of view, we 
notice that the dynamics and in particular, the structure of pe- 
riodic orbits are not affected at all by the perturbation. In this 
context, the question that immediately arises is how the per- 
turbation has to modify the dynamics in order to guarantee 
the required independence. We develop the answer within the 
short periodic orbit approach [20], where the eigenfunctions 
of a chaotic system are described in a scar function basis set 
lEltl : these wave functions are supported by the shortest pe- 
riodic orbits of the system, with periods up to the Ehrenfest 
time. On the other hand, each matrix element in this basis 
includes a phase depending on the difference of actions be- 
tween the corresponding periodic orbits [22]. So, to get sta- 



tistical independence between the perturbed and unperturbed 
eigenfunction sets, the perturbation has to modify the action 
of the used periodic orbits by at least h, in a more or less ran- 
dom way. Specifically, the set of numbers [AS ^ / {2%ti)\ (mod 
1), where pi labels periodic orbits with period shorter than the 
Ehrenfest time, has to be distributed uniformly in the interval 
[0, 1); AS^ is the action of the unperturbed periodic orbit, jj., 
minus the corresponding action in the perturbed case. 

In conclusion, our results demonstrate that quantum sys- 
tems with classically chaotic dynamics react in a universal 
way as a consequence of perturbations of classical nature. 
Specifically, the LDOS is a Breit-Wigner distribution, even 
for strong perturbations. Moreover, we derive a semiclassi- 
cal expression for its width that is accurate for paradigmatic 
systems of quantum chaos as the cat maps and the stadium 
billiard. As a final remark, we would like to notice that our 
semiclassical result reproduces in part the old one obtained 
by Wigner ^ within the random matrix theory [23]. This 
fact implies that the connection between chaotic systems and 
random matrix theory, uncovered by the cellebrated Bohigas- 
Giannoni-Schmit conjecture 12411 . is stronger than believed. 
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